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Abstract 

The QCD light quark mass renormalized at a 1 GeV scale in the 
MS scheme is obtained from the numerical results of the lattice QCD 
simulation with staggered fermions. The primary emphasis is given to 
the connection between the lattice and continuum parameters. The 
results are compared with those from the QCD sum rule. 



1 Introduction 



A connection between lattice and continuum renormalized parameters is 
necessary to compare the lattice observables with the continuum observables 
(or experimental data) and to check whether the numerical simulation of 
lattice QCD makes physical sense. The bridge for the QCD coupling between 
the lattice and the continuum was done through the weak coupling expansion 
to one-loop order[l, 2, 3, 4]. The bridge for the QCD quark mass with 4-flavor 
staggered dynamical fermions was done in Ref. [5, 6]. 

Building on that work, we will attempt to relate the bare, light quark mass 
of 1.6 MeV obtained from 2-flavor staggered fermion simulations [9, 8, 7] with 
the renormalized light quark mass of between 3.5 and 9 MeV defined at a 1 
GeV scale in the MS scheme and deduced from the QCD sum rules[10]. 

In order to make this comparison we must assume that the coupling 
constant used in the lattice calculation (/3 = 6/g 2 = 5.7) is within the per- 
turbative region[ll, 12]. 

We use and compare a number of possible perturbative approaches in- 
cluding the mean field method (tadpole improvement) suggested by Lepage 
and Mackenzie. 

2 Lattice QCD with Staggered Fermions 

Lattice QCD has a lot of difficulties in implementing quark flavor dy- 
namics. There are two popular methods to put the fermions on the lattice: 
one is the Wilson fermion formalism and the other the staggered fermion 
formalism[5]. In the limit of zero quark mass, the staggered fermion scheme 
has remnants of chiral symmetry and may be preferred for the numerical 
calculation of the meson and hadron spectrum. 

2.1 Lattice QCD Action 

The current numerical simulation of lattice QCD is based on the following 
action [7, 13, 14]. 

X{x)ul{x)x{x + a„) - x(x + a^U^xix) 
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-maJ2x( x )x( x ) + Sgi uon 

X 

= -X*(p+ma)X + S g i uon 

where 

and the definition of r]^(x) is given in Ref. [5, 15]. The current approach to 
flavor dynamics is as follows: 

^ = fd[U] exp(-S gtllon )[det(|) + ma)]% ^ 



Nj_ \~J 

J d[U] exp(-S g i uon ) [det (p+ ma)] i 

for Nf degenerate flavors of mass m. It is important to notice that m is the 
bare quark mass and an input parametric mass for the numerical simulation. 
We will let m = M (a) represent the physical value for this bare lattice quark 
mass for a given value of the lattice spacing a. 

2.2 Bare Quark Mass on the Lattice 

From the numerical results [8], at (3 = 5.7, 

Mia 2 = -0.011(3) + 6.48(14)ma (% 2 = 4.2) 

M p a = 0.31(2) + 9. 7(8)ma ( x 2 = 0.06) 

M N a = 0.47(3) + 15.0(1.2)ma (x 2 = 2.6) 

Ua = 0.044(1) (3) 

We can choose our conventions in the following way. The lattice QCD scale 
is obtained such that M p has the physical value in the limit of m — > 0. The 
bare quark mass, M (a) is obtained such that the ratio, ^ has the physical 
value in the limit of m — > M (a). At (3 = 5.7, 

- = ^ = 2ASGeV 
a 0.31 

0^Ml = 

6.48 aM 2 p y 1 
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Since Mfj/Mp = 1.52 in the calculation of Ref. [8], we will get different values 
of 1/a and M (a) if we choose Mjy to set the scale; 

- = ^- = 2MGeV 
a 0.47 

47 2 M 2 

We may also choose f w to set the scale 1/a. 

- = -^ = 2.12GeV (6) 
a 0.044 ^ ; 

3 Mass Renormalization in the Continuum 



The most common renormlization scheme in continuum QCD is the MS 
scheme in Feynmann gauge. Using dimensional regularization, the one-loop 
self-energy contribution can be calculated as follows: 



d d k „ 1 Sj, 

i(j> + fi) + M v P 
= 2^E 1 + M S 2 (7) 
and d = 4 — e , 



where M and g are the bare quark mass and coupling respectively. The 
renormalized couplings are related as follows: 



e/2 

go == 



Z rns = 1 _ 9rns^ 2 fll _ 1 

9 16vr 2 e V 2 3 * 



The quantities in Eq.(7) are given by 

N 2 - 1 
R(N C ) = " 



2N C 

-.2 



Si = 9 f^R{N c ){ 1 -- f d x {l- x )H^]- 1 -) (8) 
8?r \e Jo 2 1 
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s 2 

A 



f±R(N c )(--2 /^xlnfA] 
x(l - x)p 2 + xMq . 



(9) 
(10) 



The renormalized quark propagator, G is related to the bare quark propaga- 
tor, Go as follows: 

G {jPi M. ms , Cjms , l^ms) 

= Z™ s GQ 1 {p,M ,g ,e)\ — ./j^ (11) 

t- iiao — ^ rn lvl ms > i/0 — Harris ^ g ifrrts 

where Gq M , g , e) = 2^ + M + M , g , e) (12) 



In the MS scheme[18], the renormalization constants Z™ s and Z™ s are 
given to the lowest order in g^j by 



ryms 



<?— 1 

^R(N C )- 
8vr 2 ^ ; e 



8vr 2 



so the renormalized quark propagator for N c = 3 is 



C 1 (p, cfaj, fi^j) = 1 



Q— 

orris 

6vr 2 



/^(l - x)ln( 
Jo 



l + a .)ln(_A 



A . 



/4u 2 



(13) 



In perturbation theory, we interpret the location of the pole in the renor- 
malized quark-propagator as the physical mass of the particle, regardless 
of the renormalization scheme[20]. So if we define the pole location as 
p 2 = —M 2 hy , Mphy should satisfy the following condition. 



M, 



phy 



l + x)ln(-=- 



4 



(14) 



-M 



phy 



Since the pole location is independent of the renormalization scheme, it will 
be important for connecting the continuum and lattice schemes below. 
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4 Mass Renormalization on the Lattice 



Since the ultra-violet divergences are already regularized by the lattice, 
only a subtraction prescription is needed for renormalization. One of the 
most important characteristics common to both MS and MS schemes is 
that the renormalization constants are independent of any dimensionful la- 
grangian parameters such as the mass[19]. We will follow a similar procedure 
in defining a renormalization scheme on the lattice, choosing renormalization 
constants in the lattice-regularization formalism that are independent of any 
dimensionful lagrangian parameters such as the mass. 

4.1 Renormalization Prescription 

In order to implement the minimal subtraction idea in the lattice regu- 
larized theory [21, 6], only divergent terms such as (lnfa/z^] — Constant) 1 with 
i < n will be subtracted consistently in the lattice regularized Feynmann 
diagrams of g^ 1 order. The arbitrary constant will be chosen in a physically 
reasonable way, as is described in detail later. In such a minimally-subtracted 
theory, it is well known that the renormalization group functions, /3(g) and 
7m(s0 are independent of the covariant gauge choice[16]. 

Let the renormalization constants, Zf 1 , Z£ and Z^f, be defined as follows 

i g i Tin. ip 



where t = — 21n(a/Xi) and is the lattice renormalization scale introduced 
by our subtraction. 

For the coupling constant renormalization, we choose a prescription so 
that Ai, given by the conventional formula [1, 2, 4] 



can be expressed by the identical formula when go(a) is replaced by gL^L,)- 



[6]: 




(15) 
(16) 
(17) 
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This requires that our choice of Z^ satisfy: 

Z^(t = 0,g L ^ L = -)) = l . (18) 

Eq.(18) insures that the renormalized coupling, Ql{^l) and the bare coupling, 
go(a) agree for = 1/a A^,. By choosing the prescription as above, we 
have defined in the conventional way [1, 2, 4]. 

Having made a choice for as above, we know that the ratio, Hl/ ^ws 
is equal to Al/A^j in order to keep the coupling common in both renor- 
malization schemes. In this sense, a renormalization scale on the lattice, 
jiL = 1/a ~ 2 GeV (Eq.(4), Eq.(5) and Eq.(6)) corresponds to about 
l^ws ~ 90 GeV [23] for the /3 = 5.7 numerical simulations [7, 8, 9]. 

The natural question is what choice of is best to do perturbation. 
Experiment suggests that the perturbative expansion converges most rapidly 
if we choose the renormalization scale equal to the energy-momentum of 
the actual physical process in the MS scheme rather than in other schemes 
{i^tki = [20])- Hence we know that fi^j = physical energy-momentum 
in the MS scheme is the best to do perturbation [20]. When working with 
this quark mass parameter, M (a) initially defined at the lattice scale a, it 
is reasonable to use fi^j = 1/a to set the physical energy scale. S 
GeV is chosen, which corresponds to % 50 MeV on the lattice. The point 
is that it is better to choose fi^j = 1/a than = 1/a in order to improve 
the convergence of the perturbation series. 

For the mass and wave function renormalization constants, we will leave 
our prescription quite flexible requiring only: 

Z^t = C 1} g L ) = l , Z$(t = C 2 ,g L ) = l (19) 

where the uncertainties, now transferred to C\ and C 2 will be determined 
later in Section 5. 

4.2 One-Loop Self-Energy and Mass Renormalization 

The Feynmann rules and gauge-fixing(Feynman gauge) we use are con- 
sistent with Ref. [4, 5, 15]. Only the two Feynmann diagrams in Figure 1 can 
contribute to the one-loop self-energy, as follows: 

£(p, M (a),g ,a) 
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+- d 4 k \{ 1 + cos(2ap M + k^) 



£„ ilv sm(ap v + k v ) + aMo^l £ a sin 2 (^ 



"| E *7p sin(a Pp ) £* ^ 4 ^ 



J 2 



(20) 



Here 7V C is the number of colors and R(N C ) is given in Section 4. The 7^ are 
the normal 4x4 Euclidean Dirac matrices and the expression in Eq.(20) is 
independent of color and flavor. It is the same for Nf = 4 as the propagator 
derived in Ref.[5]. and £3 in Eq.(20) are 

= ^R(Nc) ( 8vr 2 r - ln(a/i L ) - /' dx(l - x) In 4 ) + °(°) ( 21 ) 

87r V 70 L^lJ/ 



and t = -0.044566, a = 0.19745 [22]. 



f22l 



In the above one-loop calculation, it is assumed that 1/a ^> p. As in 
Eq.(ll) and Eq.(12), the renormalized propagator, Gl is related to the bare 
propagator, Go, as follows: 

G^^Ml^l,^) = Z^G- 1 {p,M ,g ,a) \m =z^m l , go =z$ gL ( 23 ) 
where 

Go\p,M (a),g ,a) = £ 27^ + M (a) + M (a), g , a) (24) 

From the presciptions in section 4.1, we can determine the renormalization 
constants up to g\ order, as follows: 

Zi = l + ^c)(%L)ly) 

ft = l + ^ji2(JVc)(31n(a ML ) + ^) , (25) 
so that the renormalized quark propagator for N c = 3 is 



1 + 



6vr 2 



8vr 2 r + C 2 



/ 1 rfx( 
Jo 



1 - aOln 



7 L 



7 



it* + M L jl + ^ (^(a - r) + ^± - jf* + x) In 



A 



If we suppose that the pole is located at p 2 = —M 2 hy , M p h y should satisfy 
the following condition: 



M phy = M L \l + -^- ( 8ir 2 (a - r) + ^ - / dx(l + x)ln 



6vr 2 v 

.2 71^2 



f 1 dx(] 
Jo 



A_ 

7 L 



(26) 



where A = A(p 2 = -M 2 hy ) . 



Eq.(26) will be used later to relate Ml and M^j defined in the two different 
schemes. 



5 Connection between the Lattice and Con- 
tinuum Parameters 

The bridge condition will be constructed, exploiting the similar form the 
renormalization group equations take in the two schemes through one-loop 
order and choosing the same pole structure for the quark parpagator in both 
schemes. 

5.1 The Bridge Condition 

In the framework of a mass-independent renormalization scheme, the 
renormalization group determines the renormalized coupling as a function 
of only /j,/ A and the renormalized mass as the product of the renormaliza- 
tion group invariant mass, M and a function of only /x/A [10], where A is an 
integration constant chosen by standard convention. Then the following two 
conditions can be chosen: 

• Coupling relation : The renormalized couplings are the same for the 
introduced renormalization scales. For example, 

5l(^l) = gmsififna) = 9 <^ T~~ = T~~" (27) 

because using the conventional definition of the A parameter in both 
schemes we have 



5l(^l) = /(t^) and fer = f{^] 
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for the same function /() through two loops. 



• Mass relation : The renormalization group invariant mass is common 
to any renormalization scheme that is independent of dimensionful la- 
grangian parameters. For example, 

^ = ^ and common M <j=* M L (/z L , M) = M™(/w, M) (28) 
again because 

M L (fi L ,M) = Mh(^) and M^(^,M) = Mh(^) , 

through one-loop order. 

In particular, for (3(g) = — flog 3 — (3\g 5 — ^g 7 — • • •, it is a well-known fact 
that the first two coefficient s(/3o and (3\) are independent of the renormal- 
ization scheme while all the other coefficients depend on the renormalization 
scheme. However for "f m {g) = "fog 2 + 7i5 4 + ' ' 'j om y the first coefficient 70 
is independent of the renormalization scheme[20]. In the framework of the 
bridge condition, only the scheme-independent parts of the (3 and "f m func- 
tions are considered and the scheme-dependent parts ignored so that /(^) 
and h(^) are universal i.e. independent of the renormalization scheme up to 
the order of our present calculation. The constant C\ introduced in Eq.(19) 
will be determined only up to one-loop order. ■^ L - is estimated up to one-loop 
order[l, 2, 4] and has no higher order corrections [3]. 

The location of the pole in the two-point Green's function should be 
independent of the renormalization scheme[17] and determines the physical 
mass of the particle. The location of the pole in the quark propagator is 
independent of the renormalization scheme. The constant C\ can be chosen 
explicitly to conform to the conventions of Eq.(28), by demanding that the 
pole location in the quark propagator is the same for both MS and lattice 
renormalization scheme. 

Therefore from Eq.(14), Eq.(26), Eq.(27) and Eq.(28), 



M L (^M) g^ 1 3^ 

M^i^M) 6vr 2 V 1 ; 2 2 

so d = (t - a) h 2 In 



JlL_ 
f^ms 



1 



3 LA 



(29) 
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where t and a are given in Eq.(22). Various values of C\ for Nf = 0, 2, 3, 4, 
5, 6 and 8, are given in Table 1. 

5.2 Renormalized Coupling Constant 

There are a number of ways to obtain the renormalized coupling con- 
stant, g(fJ-). Once we know the QCD A parameter, we know the coupling 
constant g(fJ-). For example, one may use the experimental value of the A 
parameter in the MS scheme [31]. For a given renormalization scale /x^jtt, 
we can obtain the renormalized coupling constant, Ql{^l) since the ratio, 
Ki,/ Kj^g is well-known [1, 2, 3, 4]. Another method uses the bare coupling 
on the lattice in order to obtain the renormalized coupling. In this article, 
the latter method is chosen in order to do everything consistently in terms 
of the lattice parameters and the lattice observables. 

From Eq.(15), the renormalization group equation can be derived, as 
follows: 



Eq.(30) can be solved, with Eq.(18) as a boundary condition [6], up to 




(30) 



where /3(g) 



/V - Pig 5 




the leading order in g 2 L , 




(31) 



where 




(32) 



From Eq.(31), for various Nf and fi^j, g^j can be obtained. 
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5.3 Renormalized Mass 



We have already obtained an expression for the mass renormalization 
constant, in Eq.(25). Thus the naive relationship between the renormal- 
ized quark mass in MS and the bare quark mass in the lattice formulation 
follows from Eq.(16), Eq.(25), Eq.(27), Eq.(28) and Eq.(29): 



-^ms( [l"ms ) 
-^ms( [l"ms ) 



9 -j^ (l^afi^) + In ^ 



As described in section 4.1, it is reasonable to use /x^ 
for the optimal perturbative expansion. From Eq.(33), 



Mrns{l/a) 



1 



8vr 2 



2vr 2 



t - a 



V 



M (a) 



M (a) 

M {a) (33) 
1/a as condition 

(34) 



Solving the renormalization group equation, we can obtain the renormalized 
quark mass at any other scale in the following way. 



lm{g) = 
where 70 



log 2 + 715 4 + •••• 



(35) 



2vr 2 



We can solve Eq.(35) such that the renormalization scale runs from 1/a to 
an arbitrary scale, Then the solution to Eq.(35) is 



TO 
2/9 



Mrnsil/a) 

„2 / o^2 



2vr 2 



t - a] 



M„(a) (36) 



Other authors quote values for M^j(/j,^j) evaluated at fi^j = 1 GeV. From 
Eq.(36), we can obtain the renormalized quark mass at that 1 GeV scale in 
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the MS scheme(Table 1), although we are entering a region where perturba- 
tion is unlikely to be accurate. 

Equation (33) is just the lowest order relationship and can be improved 
through leading logarithmic approximation. From Eq.(16), the renormaliza- 
tion group equation is derived as follows: 

(-2^ + /3(<7l)^ -7™Gfc)) Zi{t,g L ) = (37) 

The solution of Eq.(34) with Eq.(19) as a boundary condition is, up to the 
leading order in g 2 L [6], 

Z^t,g L )= [l+A^Ct-d)] '*> 
and from Eq.(16), Eq.(27) and Eq.(28), 

1 TO 

M^- A \^) = j^Moia) = [1 + (3 gUt ~ CiFM (a) (38) 

m 

where the renormalized coupling and t can be obtained from Eq.(31) and 
Eq.(32). As one can see in the above, up to g^j order, Eq.(38) is identical 
to Eq.(33). From Eq.(36) and Eq.(38), for various Nf, for (3 = 5.7 and for 
firns = 1 GeV, M^s and M^- A - are given in Table 1. 

5.4 Renormalization Group Invariant Mass 

Using the conventions of Ref. [10] and expanding the function, 7771(5), the 
renormalization group invariant mass is defined, up to leading order in g 2 L , 

by 

Mr 

M = V • (39) 

From Eq.(15), Eq.(16), Eq.(31), Eq.(38) and Eq.(39), M is related to M (a), 
as follows: 

" = \Mw " Cl ) "° M ° (a) (40) 
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Eq.(40) is independent of the renormalization scale, or fi^j, which 
means that M is less sensitive to the renormalization scheme. It is reasonable 
to choose the renormalization group invariant mass as the physical quantity 
to extract from a lattice calculation, since we don't need to introduce the 
renormalization scale, or The renormalization group invariant quark 
mass is given in Table 1, from Eq.(40). 



5.5 Light Quark Mass for N f = 2, (3 = 5.7 

In this section, it is explicitly explained how the renormalized mass at the 
given scale (1 GeV in MS scheme) and the renormalization group invariant 
mass[10] can be estimated in the manner described above using the numerical 
results of the lattice QCD simulation. 

Let us start with the renormalized coupling, g^rj(/J,^rj) or Ql{^l)- We know 
that (3 = 5.7 corresponds to the bare coupling constant, g^(a) = 1.05263 
since gg(a) = 6//3. From Ref.[4] and Table 1, it is known that Al/A^j = 
0.02279 for Nf = 2. Since = (Al/ A^j)/j,^j, we know that fi^j = 1 GeV 
corresponds to \i L = 0.02279 GeV. From Eq.(4), Eq.(5), Eq(6) and Ref.[7], 
a' 1 = 2.23 GeV and t = -21n(a/z L ) = 9.1707. For N f = 2, /3 = ^ff,, 
which is defined in Eq.(30). Since we know go(a), t and /3 , the renormalized 
coupling constant at a scale of 1 GeV is from Eq.(31), 



= lGeV) = gl(p L = 23MeV) = 2.57 . (41) 

In a similar way, we can obtain the renormalized coupling constant at a scale 
of 1 1 a: 

sLC/x™ = 1/ a) = g^ = 51 MeV) = 2.05 . (42) 

Now let us go ahead to obtain the renormalized light quark mass at a 
scale of I /a. We know g^g{l /a) and (t — a) from Eq.(22). The bare quark 
mass M (a) = 1.6 MeV at most from Eq.(4), Eq.(5) and Ref.[7, 9]. So the 
renormalized light quark mass at a scale of 1/a is naively, from Eq.(34), 

M (a) 

M^(l/a) = M L (51 MeV) = 2.69MeV. (43) 
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We can also obtain the renormalized quark mass at a scale of 1 GeV from 
Eq.(41), Eq.(42) and Eq.(36). 

M^(l GeV) 

[\ , , 1 = 1-84 
M (a) 

M^(l GeV) = M L (51 MeV) = 2.95 MeV . (44) 

The renormalized light quark mass at a scale of 1 GeV, improved by 
leading logarithmic approximation is, from Eq.(38), 

M^- A \\ GeV) = Mf- L - A (51 MeV) = 1.64 M (a) = 2.62 MeV . (45) 

Now let us obtain the renormalization group invariant mass. Knowing 
Al/A™ = 0.02279 from Table 1 and Ref.[4], we can determine d = -5.5098 
from Eq.(29). Since we know 70, /3o, go(a) and C\, the renormalization group 
invariant mass is, from Eq.(40), 

M = 3.52M (a) = 5.63 MeV . (46) 

Eq.(33) is just a crude result of one-loop renormalization, while Eq.(40) 
is more reliable since it includes all the leading logarithmic contributions and 
is independent of all the renormalization scales introduced, which means less 
dependence on the conventions. The 10% difference between the numbers in 
Eq.(44) and Eq.(45) suggests the size of the omitted higher order corrections. 

6 Mean Field Theory 

Lepage and Mackenzie showed that tadpole diagrams are the main source of 
the large difference between the bare lattice coupling, go(a) and the renor- 
malized coupling, ffMsifJ-'Ms = [1^]' T ne y su gg es t a meanfield method for 
removing the dominant effect of tadpole diagrams. Here we apply their mean- 
field method to the staggered fermion formalism, in order to improve the 
estimation of the renormalized quark mass on the lattice non-perturbatively. 

6.1 Tadpole Improvement for the Coupling Constant 

The matching of lattice operators with continuum operators is based on the 
expansion 

U^x) = e ia9A ^ -> 1 + zagA^x) (47) 
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when the lattice spacing a is small. But higher order terms in the expansion 
of U,j, contain addtional factors of agA^ and the contraction of A M (x)'s with 
each other generates ultra-violet divergence (i.e. oc ^) which cancel the 
addional powers of a (i.e. oc a n ). These terms are of the order of g n , not 
suppressed by the powers of lattice spacing a, and are very large. These are 
called the QCD tadpole contributions[12]. 

These large tadpole contribution causes poor perturbative expansion in 
lattice QCD. Lepage and Mackenzie suggested the mean field method[12] in 
order to refine naive perturbative expansion by removing tadpole contribu- 
tion. They noticed that the vacuum expectation value of the link matrix is 
smaller than 1. They suggested that the appropriate connection with the 
continuum gauge field is 



u (l + iagA^x)) 



(48) 



where u represents the mean value of the link. They choose u in a gauge 
invariant way: 



u 



Re < \ttU d > 



(49) 



Other choices of u are possible and equally good. It is the rescaled links 
that should be expanded around unity. The staggered fermion action given 
in Eq.(l) can be rewritten in terms of ip = -v/uoX 



' gluon 



\~nvi x 

M 



u 



X iP(x + <v) - i>(x + aj^^-^(x)} 



Uq 



Uo 



a^2ip(x)ip(x) + 



E/3 



MF 



1 ReTrffp 

%. ut 



(50) 
(51) 



where 



/3. 



MF 



2N C 

2 ' 

9mf 



2 

9mf 



9 2 o{a) 



Re < |TrC/ D > 



The perturbative expansion for Trt/ n given in Ref. [32] is 



Re < -TrUn > = 
3 



1 - \gl{a) + 0{gt{a)) 



(52) 



(53) 
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We can relate g 2 MF to g 2 perturbatively as follows: 



q2mf = Re 4 ( Trk > = + ^° (a) + °^ a ^ 

From Eq.(31), Eq.(32) and Eq.(54) we can also relate g 2 MF to g^j(fJ-7fn) per- 
turbatively: 



~\g 2 MF + °{9mf)] (55) 
for Nf = 2, s£u(/Xmj) = 5mf[! ~ /3o5mf{ 2 m ( a M^) ~ 2.117} 

+0(g 4 MF )] (56) 

W™ = f, g 2 -C-) = g 2 MF [l - 0M06g 2 MF + 0(g 4 MF )] (57) 

As one can see in Eq.(57), g 2 MF is extremely close to g^j{^), where the 
difference between these two coupling constants is only a few percent from 
the standpoint of pertubative expansion. This means that by running the 
renormalized coupling in the MS scheme at a scale of ^, one can absorb the 
large tadpole contribution into the renormalized coupling constant. 



6.2 Tadpole Improvement for the Quark Mass 

As the lattice spacing, a goes to zero, the action given in Eq.(50) will be 

f - M 

S = - / V(z)[£>-7 + — + Sglvon (58) 

Jx U 

Here the claim is that ip is better matched to the continuum quark field than 
X- From the above action, we notice that Mmf = takes the place of 

the bare quark mass on the lattice. 

Hence the renormalized quark mass can be obtained as follows: 

M^fi^j) = Z m (afi 7fri ,g^ ri (fi 7 f r i))M (a) 
= Z m u M M F(a) 

= Z m {{aii^,g 2 MF )MMF{a) (59) 
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where Z m = Z m u . 

From Eq.(33), Eq.(49), Eq.(53), Eq.(54) and Eq.(55), the perturbative 
expansion of Z m can be obtained as follows: 

■^m(( a Mms) 9mf) 



1 



r 



12 



2 

5mf 



2vr 2 \ 



ln(a/x^) + 



8vr 2 



T - cr 



1 vr^ 
6 + 6 , 



(60) 



The numerical data for u is obtained by the Columbia lattice group at 
= 5.7 for N f = 2 [9]: 



Re < ]-TrUn >= 0.58 



(61) 

u = 0.87 . (62) 

Putting everything altogether, we can obtain the renormalized quark mass: 

(63) 
(64) 

(65) 



M M f{o) = IMMeV 
Z m (l,g 2 MF = 1.815) = 1.45 

AW-) = 2.66MeV 



This final result (Eq.(65)) is in good agreement with that of Eq.(43) in the 
section 5.5. This good agreement means that in fact through perturbative ex- 
pansion with respect to <fe?(l/a), one can remove the tadpole contributions, 
which is very large when one does the pertubative expansion with respect to 
the lattice bare coupling, </q. 



7 Light Quark Mass from QCD sum rule 

There are various methods to determine the light quark mass in contin- 
uum QCD[10, 25, 26, 27]. The actual determination of the light quark mass 
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in the framework of QCD can be done reasonably through the QCD sum 
rule formalism by Shifman, Voloshin and Zakharov (S.V.Z.) [28, 29]. In that 
formalism, they study the two point correlation function of the divergence of 
the axial current (isospin sector). They determine the running quark mass 
at 1.2 GeV scale with Nf = 3 i.e. three dynamical quarks moving around. In 
the MS scheme, they choose renormalization scale fi^j to be the same as the 
energy-momentum scale of the physical process. Their choice, = 1 ~ 1.2 
GeV means physically that the charm and bottom quarks are decoupled [30] 
and that only three light quarks contribute to the dynamics of the QCD. In 
other words, at fi^j = 1 ~ 1.2 GeV, we can not decouple the strange quark 
from the dynamics of the QCD. Now we have a problem that lattice QCD 
is simulated numerically with 2 dynamical quarks but the QCD sum rule 
assumes 3 dynamical quarks. 

From the recent experimental data[31], we have the following QCD A 
parameter: 

AM N f = 3) = 250 ± 150MeV 

where the error is partly a combination of statistial and sytematic errors and 
partly due to the scale uncertainty. We should use A^(Nf = 3) at the scale 
below the charm quark mass. Using A^j(Nf = 3) = 250 ± 150 MeV, the 
QCD sum rule gives the following light quark mass[10]. For Nf = 3, 

M^{^=lGeV) = Mu + Md = 3.5 ~ 9.0 MeV 
M = 4.0 ~ 13.0 MeV , 

where the larger value corresponds to the smaller < 150 MeV and the 
smaller value corresponds to the larger A^ > 300 MeV. 

The QCD sum rules take into consideration the resonance contribution 
but not the continuum contribution to the imaginary part of the two point 
correlation function of the axial current divergence. Omitting this contin- 
uum contribution from the non-perturbative low-energy region can introduce 
around 20% errors to the above expectation value of the light quark mass[28]. 

Even though there is a difference in the number of dynamical quarks, 
the lattice QCD expectation value of the quark mass with Nf = 2 given in 
Section 6.5 agrees with that of the S.V.Z. QCD sum rule with Nf = 3 in the 
above. 
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But in order to do the exact comparison, the lattice QCD simulation with 
three dynmical quarks (three sea quarks: one of them is 10 ~ 20 times heavier 
than the other two light sea quarks) is necessary since the strange quark can 
not be decoupled from the QCD dynamics at the energy-momentum scale of 
1 GeV[30]. 

8 Conclusion 

The QCD dynamical quark mass is renormalized in two different renor- 
malization schemes (MS and lattice-regularized minimal subtraction). The 
bridge conditions are chosen to make connection between the two schemes. 
The ratios of the continuum renormalized quark mass to the lattice bare 
quark mass are given in Table 1. 

Ref. [7] contains an earlier attempt to compute the renormalized quark 
mass. Their value of 2.2 MeV is somewhat different from that of this arti- 
cle(2.95 MeV). Fukugita et al. use the one-loop renormalization and Eq.(3.7) 
in Ref. [7] is the same as Eq.(33) except for the coupling constant that ap- 
pears. This 25% difference comes from the different choice of the coupling: 
the bare coupling on the lattice is used in Ref. [7] while the renormalized 
coupling is used in this article. It has been pointed out that the lattice bare 
coupling constant may be a poor choice as an expansion parameter and that 
the use of the improved coupling constant including renormalization due to 
gluon tadpole contributions improves the reliability of the lattice perturba- 
tive expansion[12]. Using the mean field method suggested by Lepage and 
Mackenzie, the tadpole-improved renormalized quark mass is obtained, which 
is extremely close to that obtained by the use of the renormalized coupling 
in the MS scheme at ^ scale. It has been proposed in other contexts [20] to 
use the renormalized coupling (in the MS scheme with the renormalization 
scale equal to the physical energy-momentum) instead of the bare coupling. 
At the least, the differences between the results obtained here and those of 
Fukugita et al. represent the size of the perturbative errors. 

The QCD sum rule predicts that the renormalized quark mass (at a 1 
GeV scale in the MS scheme) is 3.5 ~ 9.0 MeV and that the renormalization 
group invariant mass is 4.0 ~ 13.0 MeV for Nf = 3, with 20% uncertainty 
[10, 28, 29]. The renormalized quark mass (at a 1 GeV scale in the MS 
scheme) obtained from the lattice QCD simulation is 2.95 MeV and the 
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renormalization group invariant mass 5.63 MeV for Nf = 2. 

The difference between 2-flavor dynamics and 3-flavor dynamics with one 
of the three flavors much heavier may be presumed to be so small that the 
comparison between 2-flavor and 3-flavor dynamics may make physical sense. 

There are two other sources of error in the lattice QCD simulation: one 
is finite-temperature effect and the other finite- volume effect. In the hadron 
mass calculations on 16 3 X 16, 16 3 X 32 and 32 3 X 32 shows approximately 
5 ~ 10% effect [33]. These small effects are completely neglected in this 
article. But the systematic analysis of finite size effects should be done to 
look into physics on the lattice more precisely. 
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A L 


Ci 


M^l/a) 

M 


M^(lGeV) 

M 


MlL± A (lGeV) 

777.S ^ > 

Mo 


M 

Mo 





0.0347 


-6.35 


1.69 


1.86 


1.62 


2.97 


2 


0.0228 


-5.51 


1.68 


1.84 


1.64 


3.52 


3 


0.0176 


-5.00 


1.68 


1.84 


1.65 


3.92 


4 


0.0131 


-4.40 


1.67 


1.83 


1.67 


4.45 


5 


0.00922 


-3.70 


1.67 


1.83 


1.68 


5.17 


6 


0.00608 


-2.87 


1.67 


1.82 


1.70 


6.19 


8 


0.00197 


-0.61 


1.66 


1.81 


1.74 


10.26 



Table 1: Here we use (3 = 5.7, a 1 = 2.23GeV . The ratio, A^/A^ is obtained 
from Ref.[4]. 
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Figure 1: One-Loop Contribution to the Self- Energy (a) gluon exchage dia- 
gram (b) gluon bubble diagram 
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